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Abstract 

Production of the K*(892)° and 7T(892)° resonances was studied via their K + ir and K ir + 
decay modes in central Pb+Pb, Si+Si, C+C and inelastic p+p collisions at 158^4 GeV (y / ijv r /v 
= 17.3 GeV) with the NA49 detector at the CERN SPS. Transverse momentum and rapidity 
distributions were measured and total yields were estimated. The yield of K* exceeds that of 
K* by about a factor of two in nucleus-nucleus reactions. The total yield ratios (K*) / (K + ) and 
(K ) I (K~) are strongly suppressed in central Pb+Pb compared to p+p, C+C and Si+Si collisions 
in agreement with the expected attenuation of these short-lived resonance states in the hadronic 
phase of the expanding fireball. The UrQMD model, although incorporating such a scenario, does 
not provide a quantitative description of the experimental results. The statistical hadron gas 
model assuming the same freeze-out parameters for stable hadrons and resonances overestimates 
the {K*)/{K) ratios in central Pb+Pb collisions by about a factor of 2.5. 

PACS numbers: 13.85.Ni,25.75.Dw 
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I. INTRODUCTION 



High-energy collisions of heavy nuclei produce a transient state of extreme energy and 
matter density in which quarks and gluons are probably briefly deconfined jl-3]. Production 
of entropy and of s, s quarks is believed to occur at the early stage of the collision and this 
process is expected to be sensitive to the phase of the created matter 4], |5] . The high-density 
state evolves into a hadron-resonance gas which finally decouples into the observed hadrons. 
The A'* (892) and K (892) resonance states contain an s and s valence quark, respectively, 
and are therefore sensitive to the level of strangeness production. However, resonance states 
have lifetimes similar to that of the fireball and may interact in the dense medium in which 
they are produced. Their mass and width could be affected {f| and scattering processes 
might destroy or regenerate them. Furthermore, daughters of those K* that decay inside 
the fireball may rescatter resulting in a changed invariant-mass spectrum. Thus the yields 
contained in the K* mass peak were conjectured to be sensitive to the duration and properties 
of the hadronic fireball stage 

Studies of K* (892) production at mid-rapidity in Au+Au, Cu+Cu and p+p collisions 
at RHIC energies were pe—bytheSTAReoHaWtionBa. This paper reports 
measurements of K* (892)° and K (892)° resonance production via their K + n and K 7r + 
decay modes at the CERN SPS in central Pb+Pb, Si+Si, C+C and inelastic p+p collisions at 
158A GeV (^snn = 17.3 GeV). Preliminary results were presented in [10[. Section II briefly 
describes the NA49 detector. Section III discusses the analysis procedure. Distributions of 
transverse momentum and center-of-mass rapidity y as well as total yields are presented 
in Section IV. These results are compared to predictions of the ultrarelativistic quantum 
molecular dynamics (UrQMD) model ll| and a statistical hadron gas model (HGM) 12j in 
Section V. The paper ends with the summary Section VI. 



II. DETECTOR 



The NA49 experimental apparatus 13[ consists of four large-volume time projection 
chambers (TPC). Two of these (VTPC) are placed in the fields of two super-conducting 
dipole magnets. The other two (MTPC) are positioned downstream of the magnets and 
are optimized for high-precision measurements of the ionization energy loss dE/dx with a 
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resolution of about 4%. The particle identification provided by the dE/dx measurement is 
complemented in the mid-rapidity domain by a measurement of the time-of-flight (TOF) 
with a resolution of about 60 ps in two TOF detector arrays positioned downstream of 
the MTPCs. The magnetic fields were set to about 1.5 T (upstream magnet) and 1.1 T 
(downstream magnet). With the lower momentum cut employed for dE/dx identification 
the detector acceptance covers the forward rapidity region for if* (892). 

The precise transverse position of each beam particle at the target was measured by three 
pairs of small proportional wire chambers (BPD) upstream of the target with a precision of 
better than 200 /im. Lead ions of 158A GeV impinged on a thin Pb-foil target of 337 mg/cm 2 
(approximately 1.5 % interaction probability for Pb ions) which was positioned 80 cm up- 
stream from the first VTPC For the study of C+C and Si+Si collisions a 3 mm thick C 
(2.4% interaction length) and 5 mm thick Si target (4.4%) were used, respectively. The 
incident C and Si nuclei were produced by fragmentation of a Pb beam of 158^4 GeV beam 
energy 13[ and were selected by magnetic rigidity (Z/A — 0.5) and by specific energy loss in 
the BPDs. The " C-beam" as defined by the online trigger and offline selection was a mixture 
of ions with Z = 6 and 7 (intensity ratio 69:31); the "Si-Beam" of ions with Z = 13, 14 and 
15 (intensity ratio 35:41:24). The trigger selected the centrality of the collisions based on 
a measurement of the energy deposited by projectile spectator nucleons in a downstream 
calorimeter. 

For the study of p+p collisions the beam line was set to select secondary protons of 
158 GeV/c momentum which were produced in a Be target by the 400 GeV/c SPS proton 
beam. The secondary protons were identified by Cherenkov counters in the H2 beamline 
resulting in a contamination by pions and kaons of less then 10 -3 . Liquid hydrogen targets of 
14 cm (year 1996) and 20 cm (later years, 2.8 % interaction length) and 3 cm diameter were 
used. A scintillation counter S4 of 2 cm diameter was positioned about 5 m downstream on 
the deflected beam line between the two VTPCs. It was used in anticoincidence with the 
beam in order to select p+p interactions. For a detailed description of detector aspects for 
p+p collisions see 
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III. DATA ANALYSIS 



The analysis of Pb+Pb reactions is based on a high-statistics data run which recorded 
about 3 • 10 6 collisions. The trigger selected the 23.5% most central Pb+Pb collisions. The 
corresponding mean number of wounded nucleons N w 15J was calculated using the VENUS 
simulation code [l^] following the Glauber model approach, and found to be (N w ) = 262 
with a systematic uncertainty of ±5. More details on the procedure can be found in Ref. 16]. 

The C+C and Si+Si collision data are more limited in statistics. For both systems about 
45 • 10 3 events were recorded for the (15.3 ± 2.4)% and (12.2 ± 1.8)% most central C+C and 
Si+Si collisions, respectively. The corresponding mean numbers of wounded nucleons based 
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on VENUS simulations are 14 ± 2 and 37 + 3 

Results on p+p reactions are based on 1.125 • 10 6 (4.18 • 10 5 ) events collected with the 
20 (14) cm long liquid hydrogen target. The trigger cross section was 28.3+0.1 mb which 
contained about 86 % of the inelastic cross section and excludes most of the elastic collisions 
(about 1 mb remaining contamination) jl^ . 

Charged particle tracks were reconstructed from the charge deposited along the particle 
trajectories in the TPCs using a global tracking scheme which combines track segments that 
belong to the same physical particle but were detected in different TPCs. A vertex fit was 
then performed using the reconstructed tracks. Particle identification is based on measure- 
ments of the 50 % truncated mean of the specific energy loss dE/dx in the TPCs which 
provide up to 234 charge samples on a track. The uncertainty of the dE/dx measurement 
for a specific track depends on its visible length and the number of associated charge clus- 
ters. The average value of dE/dx is a universal function of the velocity of a charged particle 
(Bethe-Bloch curve) and thus depends on its mass for a given momentum. 



A. Pb+Pb collisions 

For Pb+Pb collisions the event vertex was determined using the tracks reconstructed in 
the TPCs. The resulting vertex distribution had widths cr(x) = 0.21 cm and a(y) = 0.15 cm 
in the coordinates transverse to the beam, and a(z) = 1.3 cm in the longitudinal direction. 
Events were accepted if the vertex z coordinate was within 1 cm of the nominal target 
position. With this requirement the background from non-target interactions was negligible. 
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Further cuts were applied at the track level. The distance of the back-extrapolated track from 
the fitted vertex position had to be below 5 cm in the (horizontal) bending plane and below 
3 cm in the vertical direction. Moreover, the number of measured points on the track had to 
exceed 25 and constitute more than 50 % of the geometrically possible maximum in order 
to eliminate split tracks. Finally, the track momentum fit was repeated including the vertex 
position resulting in a typical momentum resolution of <r{p)/p 2 ~ (0.3 — 7) • 10 -4 (GeV/c) -1 
depending on track length. 

Figure [U shows a density plot of dE/dx as a function of momentum p for accepted 
positively (a) and negatively (b) charged particles showing bands for various particle 
species. K and ir meson candidates were selected by requiring a momentum in the range 
3< p <100 GeV/c and a measured dE/dx value in a band of 2.5 and 3 standard deviations, 
respectively, around the expected mean values. Expected losses due to the cuts are small 
(< 2 %) and no correction was applied. Systematic biases of the K* yields from uncertainties 
in the fit of the Bethe-Bloch function are estimated to be below < 5 %. 




5 1 1 5 2 0.5 1 1.5 2 

log (p/(lGeV/c)) log i() (p/(lGeV/c)) 



FIG. 1: (Color online) Specific energy loss dE/dx measured in the NA49 TPCs versus momentum 
p for positively (a) and negatively (b) charged particles in central Pb+Pb collisions. Curves in (a) 
and (b) show the acceptance limits for K + and tt~ , respectively. 

Raw yields of the if* (892)° and K (892)° resonance states were extracted from the 
invariant-mass distributions calculated for K + ir~ and K~ir + pair candidates, respectively. 
First, invariant-mass 

m mv (Kir) = ^(E K + E n y-(f K + f n y (1) 
distributions were computed for all selected K + ir~ and K~ir + candidate pairs in the events. 
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Next, similarly obtained distributions from pairs taken from different events of the same mul- 
tiplicity class and normalised to the number of real pairs were subtracted in order to reduce 
the dominant contribution from combinatorial background. The resulting invariant-mass 
distributions are plotted in Fig. [2] for transverse momenta px < 2.0 GeV/c and the rapidity 
range 0.6 < y < 0.9. The peaks due to the i^*(892)° and K*(892)° resonance states are 
clearly seen above a strongly mass- dependent residual background. The mixing procedure 
preserves the inclusive single-particle phase space distributions, but destroys all correlations 
between particles. It therefore cannot fully describe the combinatorial background in real 
events which is presumably partly shaped by effects such as energy and momentum con- 
servation as well as reflections from other resonance states. The raw number of i^*(892)° 
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FIG. 2: (Color online) Invariant mass distribution of K + ir~ (a) and K~ir + (b) for < 2.0 GeV/c 
in the rapidity region 0.6 < y < 0.9 after subtraction of mixed-pair background in central Pb+Pb 
collisions. The fitted polynomial background is shown by the dashed curves, the sum of fitted 
polynomial background and signal Breit-Wigner function by the solid curves. 



and K (892)° were derived from fits to the mass distributions after the mixed-event back- 
ground had been subtracted. The fit function was chosen as a sum of a linear or a 2nd order 
polynomial background and a Breit-Wigner distribution: 

dN r 



dm* 



C 



TT iTUi. 



m ) 



(2) 



where m and T are mass and width of the K*, and C is a normalisation factor. Examples 
of such fits in the rapidity range 0.6 < y < 0.9 are shown in Fig. [2j The fits are seen to 
provide a good description of the mass distributions in the fit range 780 < mj„„ < 1100 MeV 
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and result in m = 886 ± 2 MeV and T = 49 ± 6 MeV for if* (892)° and m = 891 ± 2 MeV 
and T = 52 ± 7 MeV for K (892)°. The mass values are somewhat smaller than the world 
average of 895.94 ± 0.22 MeV 19J, and were found not to depend significantly on rapidity 
or pt as demonstrated by Fig. [3] (a,b) (example for the K* (892)° which provides the better 
statistical accuracy). Scaling of the magnetic field value by the upper limit of its systematic 
uncertainty of ± 1 % changes the fitted mass value by about ± 5 MeV. Thus the observed 
mass shift is at the limit of significance. The STAR experiment at RHIC also found a 
similarly reduced mass (see Fig. [3] (b)) but only for px below about 1 GeV/c [9J]. The fitted 
width agrees well with the world average of 48.7±0.8 MeV (the invariant mass resolution 
is about 6 MeV 20[). No significant variation with rapidity or px (see Fig. [3]c,d) was found 
in agreement with results from STAR {{J. 

Correction factors for acceptance and reconstruction efficiency were derived from Monte 
Carlo simulations. K* (892)° and K (892)° were generated with realistic distributions in 
transverse momentum and rapidity and then passed through the NA49 simulation chain 
based on GEANT 3.21 2l| and a specific TPC signal simulation software. These signals were 



then added at the TPC signal level to the raw data of real events (embedding). Finally, the 
hybrid events were reconstructed and analysed like real events. A matching step associated 
the reconstructed tracks with the originally generated tracks. Resulting invariant-mass 
spectra are plotted in Fig. H] and demonstrate that neither the mass peak position nor 
its width are affected by the measuring resolution. 

Correction factors for reconstruction inefficiencies, limited geometrical acceptance and in- 
flight decays were obtained by comparing the K* (892)° and K (892)° yields extracted from 
the reconstructed hybrid events to the generated yield. Resulting efficiencies for K* (892)° — > 
K + ir~ are shown as a function of rapidity and transverse momentum in Fig. [5] and range 
between 0.4 and 0.8. Values for K (892)° K~ ~tt + are the same within statistical errors. 
Efficiencies with embedding (full symbols) are lower than those obtained by simulating only 
K* (open symbols) but suggest only a modest effect of the high track density in the TPCs. 

Differential K* (892)° yields were obtained by fitting the invariant-mass spectra after 
subtraction of combinatorial background in bins of transverse momentum and rapidity. In 
order to ensure the stability of the fits the if* (892)° mass and width were fixed to the world 
averages. The number of K* (892)° in each bin was calculated as the integral of the Breit- 
Wigner function. Correction factors were then applied for the reconstruction efficiency and 
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FIG. 3: (Color online) Fitted mass values mo (a,b) and width F (c,d) of the i^*(892)° peak in 
the K + ir~ mass distribution versus rapidity y (a,c) and transverse momentum (b,d) in central 
Pb+Pb collisions. Dots show the fitted values with statistical error bars, and bands indicate the 
systematic uncertainties. The dotted horizontal lines indicate the world average values for mo and 
r [l(J]. For comparison results are shown from the STAR collaboration at RHIC 9] for Au+Au 
collisions at y/SNN= 62 GeV (triangles) and ^/snn= 200 GeV (squares). 



the decay branching ratio (66.7 %). 

Systematic errors were estimated by varying the identification criteria for kaons and the 
details of the fit procedure applied to the mass distributions. Increasing or decreasing the 
width of the dE/dx selection by half a standard deviation of the dE/dx measurements led to 
changes in the yields of around ± 7%. Extending or narrowing the mass range of the fit by 
50 MeV or changing from a linear to a 2nd order polynomial background affected the results 
by about ± 10%. Other sources of uncertainty, like using a mass value different from the 
world average or changing the inverse slope parameter T in the efficiency calculation, were 
much smaller. The total systematic error was estimated as half the range covered by the 
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FIG. 4: (Color online) Invariant mass spectra of simulated and embedded K*(892)° (a) and 
if* (892)° (b) calculated from reconstructed matched tracks. 
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FIG. 5: (Color online) Combined reconstruction efficiency and acceptance of K* (892)° — > K + ir~ 
as a function of rapidity (transverse momentum range < pt < 2.0 GeV/c) (a) and transverse 
momentum (rapidity range 0.43 < y < 1.78) (b) in central Pb+Pb collisions. Full and open symbols 
show results with and without embedding into real events (see text). The decay branching ratio is 
not included in the plotted efficiencies. 



results obtained when varying the dE/dx cuts and the fitting procedure as just described. 



More details of the analysis procedure for Pb+Pb collisions can be found in 



22] 



B. p+p collisions 

For p+p collisions the interaction vertex was determined from the back-extrapolated 
tracks and the trajectory of the individual beam particle in the target which was measured 
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FIG. 6: Examples of invariant-mass distributions for K + tt~ (a) and K~ir + (b) after subtraction 
of mixed-event background in inelastic p+p collisions. The curves show the fits with the sum of 
Breit-Wigner functions to describe the signals of the K* (892)° and K% (1430)° and the contributions 
from their reflections (see text). 

by the BPDs. In order to obtain a clean sample of p+p collisions, only events with a 
successfully fitted vertex differing in position by less than ±9 cm (±5.5 cm for short target) 
in z-(beam-)direction from the target center and having a radial distance of less than 1 cm 
from the beam axis were accepted, thus minimizing contributions from interactions in the 
mylar windows. From data taken with the liquid hydrogen removed from the target vessel 
("empty-target runs") the remaining fraction of background events was estimated to be 
below 1 % and therefore no correction was applied. 

Correction procedures were devised to obtain the yield per inelastic p+p collision. The 
efficiency and accuracy of vertex reconstruction using tracks back-extrapolated from the 
TPCs vary with the charged-particle multiplicity of the event. The efficiency was derived 
from the probability that a successful vertex fit is obtained and that the fitted location fell 
inside the cuts. The corrections amount to 30% for events with 3 tracks, but rapidly drop 
to 6% for events with 7 tracks. Furthermore, about 7% of the triggered events have no 
accepted tracks in the detector. Half of these can be attributed to the 1 mb contribution of 
elastic scattering events to the 28.3 mb trigger cross section, the other half are most likely 
due to singly diffractive events. In order to obtain the yield per triggered inelastic event, 
the number of events with tracks in the TPCs was scaled up by 3.5% for the per event 
normalisation. 

A further correction was applied for the 14.4 ± 1 % of inelastic events which do not 
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give rise to a trigger 



14j . Their contribution to the inclusive production cross section was 



found to gene rally depend on p?, Feynman x and the type of the produced particle under 
consideration [23[. However, for charged kaons only a weak rapidity and no significant px 
dependence was observed. For K* resonance production we assume a similar behavior and 
thus estimate that on average the observed K* yield should be scaled up by 5 ± 1 %. On 
the other hand, the number of events used for normalisation has to be increased by 14.4 % 
in order to obtain the yield per inelastic event. To account for the trigger loss we therefore 
scale down the measured i^*(892)° and J^*(892)° yields per event as defined in the previous 
paragraph by 10 ± 1 %. 

For p+p reactions invariant-mass spectra could be extracted with much lower combina- 
torial background than for A+A collisions. Selected kaons and pions were required to have 
more than 30 points per track, a momentum in the interval 4 < p < 50 GeV/c, a transverse 
momentum of pt < 1.5 GeV/c, and a measured dE/dx value within ±1.5 <JdE/dx around 
the expected dE/dx position. Pairs were entered into the invariant-mass distributions with 
the appropriate event- multiplicity dependent correction factors. The background in the 
invariant-mass spectrum was determined by mixing kaons and pions from different events. 
This distribution was normalized to the same number of entries as the real event spectrum 
and subtracted, resulting in a small undershoot around the i^*(892)° signal (see Fig. |6]). 
This undershoot structure is well described by simulations of invariant-mass distributions 
resulting from K*(892)° deca ys a nd the contributions of the i^*(892)° decay products to 



the mixed event background 



24] (see curves in Fig. [6]). No additional subtraction of re- 



maining background is necessary here. In the simulation an expected contribution from the 
-^(l^SO) (24 % of the fC(892)° yield j^) was accounted for. Its inclusion does improve 
the description of the invariant-mass distribution, but does not influence the i^*(892)° yield. 
The fitted mass of the i^*(892)° was 892 ± 5 MeV, consistent with the world average [3] 
The width was also found to agree with the world average. 

Yields of the K* (892)° and K (892)° were extracted by fitting the normalisation factor 
of the simulated invariant-mass distributions in bins of rapidity and transverse momentum. 
The simulation used the world average values for masses and widths and took into account 
the effect of geometrical acceptance and kaon decay in flight. Losses owing to reconstruction 
inefficiency are negligible in p+p collisions. Global corrections were applied for dE/dx par- 
ticle identification cuts and the decay branching fraction. Systematic errors were evaluated 
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by changing track cuts of the selected kaons and pions and amount to 8% for the integrated 
yield. A conservative systematic error of 4% is assigned for the uncertainties of the vertex 
reconstruction efficiency and the trigger loss corrections. The final systematic error is taken 
as the quadratic sum of all these contributions and amounts to 9%. More details on the 
analysis procedure for p+p collisions can be found in 26]. 
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FIG. 7: (Color online) Invariant mass distribution of K + ir~ (a,b) and K~ir + (c,d) in the rapidity 
region 0.2 < y < 1.8 and < 1-5 GeV/c after subtraction of mixed-pair background. Results are 
plotted for C+C collisions in (a,c) and for Si+Si collisions in (b,d). The fitted residual background 
is shown by the dashed lines, the sum of fitted background and signal Breit-Wigner function by 
the solid curves. 



C. C+C and Si+Si collisions 

The analysis of C+C and Si+Si collisions follows a strategy similar to that employed for 
Pb+Pb and p+p collisions, however being limited by low statistics. The available number of 
events was not sufficient to extract rapidity and transverse mass distributions, but allowed 
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only to estimate a total yield. K* (892)° and K (892)° invariant-mass spectra (see Fig. [7J 
were extracted in a wide range of rapidity 0.2 < y < 1.8 and transverse momentum px < 
1.5 GeV/ c. Kaons and pions were selected in the momentum range of 4 < p < 50 GeV/c and 
Pt < 1 GeV/c. The minimum number of points required per track was 50. Pions and kaons 
were identified by dE/dx within a band of typically ±1.5 standard deviations around the 
Bethe-Bloch value. Both cuts were varied in order to estimate the systematic uncertainties. 

A similar strategy for invariant-mass spectra, background calculation, and signal extrac- 
tion was used as for p+p collisions. However, subtraction of the mass distributions from 
mixed events did not completely remove the combinatorial background as can be seen from 
Fig. [7J The remaining background was parameterized as a straight line in the K* (892) mass 
region (dashed lines in Fig. [7J. The signal was then obtained by fitting the sum of this 
linear background and the if* (892) line shape (using the world average values for mass and 
width including its reflections as was done for p+p reactions. The results of the fits 
are shown as the solid curves in Fig. [7J 

The extracted raw yields were corrected for acceptance, the dE/dx identification cuts 
and the decay branching fraction. The mean acceptance of if* (892)° and if (892)° in the 
selected wide phase space was calculated by Monte Carlo simulations assuming a width of 
the rapidity distribution similar to that of charged kaons and inverse p^-slopes similar to 



that of the 0-meson [18| as its mass is close to that of the If* (892). A typical value of 
the mean acceptance requiring more than 100 points for both, the pion and kaon track, is 
14%. Varying the assumptions on background shape and the kinematic distributions of the 
if* (892) changes the mean acceptance by 10% only. Extracted yields have statistical errors 
of 15-20%. Varying background assumptions and selection criteria of kaons and pions, yields 
change by 15% at most. We thus assign a combined systematic and statistical error of 30%. 



More details on the analysis can be found in [26]. 



IV. RESULTS 



A. Pb+Pb collisions 



Yields of if* (892)° and if (892)° per event were extracted for the region of rapidity 
0.3 < y < 1.8 and transverse momentum < pt < 2.0 GeV/c. Efficiency-corrected results 
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are plotted as a function of rapidity y in Fig. [S] and listed in Table [H They include a 
small extrapolation in pt beyond 2.0 GeV/c based on the exponential parameterisation of 
the invariant pt distribution using the temperature parameter T fitted in the measured 
Pt range (see below). The rapidity distributions decrease with increasing y and suggest a 
maximum at mid-rapidity in view of the forward-backward symmetry of the reaction. 

TABLE I: Yields of K*(892)° and 7T(892)° per event in central Pb+Pb collisions as a function 
of rapidity y and integrated over transverse momentum. Both statistical (first) and systematic 
(second) errors are listed. 



y 


*(#*(892)°) 


$(tf*(892)°) 


0.3 - 0.6 


2.52 ± 0.22 ± 0.73 


1.52 ± 0.19 ± 0.45 


0.6 - 0.9 


2.71 ± 0.18 ± 0.47 


1.93 ± 0.13 ± 0.37 


0.9 - 1.2 


2.31 ± 0.16 ± 0.40 


1.21 ± 0.11 ± 0.31 


1.2 - 1.5 


2.08 ± 0.15 ± 0.49 


0.81 ± 0.10 ± 0.24 


1.5 - 1.8 


1.88 ± 0.13 ± 0.65 


0.45 ± 0.08 ± 0.26 


total yield (Gauss fit) 


10.3 ± 0.4 ± 2.5 


5.2 ± 0.3 ± 1.7 



Estimates of total yields were obtained by fitting Gauss functions to the rapidity distribu- 
tions centered at mid-rapidity. In order to better constrain the fit the width of the Gaussian 
was fixed at 1.2 units in rapidity. This choice was motivated by the observed systematics 
of the widths o 
UrQMD model 



rapidity distributions in Pb+Pb collisions [27J and the prediction of the 
The results are also included in Table HI 



For obtaining distributions in transverse momentum yields were extracted in four px bins 
for the rapidity range 0.43 < y < 1.78. The results for d £™ d are plotted in panels (a,c) of 
Fig. [9] and listed in Table [Til From these measurements the yield as function of transverse 



mass m T = \/p\ + m\ (where m is the K* mass) was calculated. The obtained values of 
dmrdy are sriown m panels (b,d) of Fig. [9] and are listed in Table IIIII An exponential 
function: 

rriT arriTdy 

was fitted to these measurements, where T is the inverse slope parameter and A a normalisa- 
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FIG. 8: (Color online) Rapidity distribution of A* (892)° (a) and if* (892)° (b) integrated over p T 
in central Pb+Pb collisions. Full symbols represent the measurements, open symbols were obtained 
by reflection around mid-rapidity. The bars show statistical errors, the bands indicate systematic 
errors. Solid curves show fits of Gaussian functions for estimating total yields, dashed-dotted curves 



depict predictions of the UrQMD model 



29] 



tion constant. The resulting values of T = 339 ± 9 MeV for A*(892)° and T = 329 ± 12 MeV 
or A (892)° are much larger than for kaons, but closer to that for the higher mass meson 



271. Moreover, — , dn , seems to exhibit a convex shape in the logarithmic representation 



tut diriTfly 

of Fi g. EE This behaviour may be due to the participation of the A* in the strong radial 
flow [28[. Indeed, a blast-wave calculation, also shown in Fig. |9] (b,d) by dotted curves, 
provides a good description of the K* spectra using parameters fitted to pion, kaon, proton 
and anti-proton spectra (T = 93 MeV, po = 0.91 28j). Alternatively, the convex shape 
could be due to the attenuation of the A* in the fireball which is expected to be strongest 
for low values of p?- 



B. p+p collisions 

For p+p reactions results are presented as yields per inelastic reaction. Transverse mo- 
mentum spectra (see Fig.CTOland Table PVT) were extracted in a rapidity range of 0.2 < y < 0.7. 
The range was chosen as close to mid-rapidity as the acceptance for the A* (892) allows. The 
my-spectra (see Fig. [TT] and Table IVII) show a thermal shape which can be well described 
by an exponential function (Eq. (13"1)) with an inverse slope parameter T = 166 ± 15 MeV 
for A*(892)° and 150 ± 14 MeV for A (892)°, respectively. These values are consistent with 
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FIG. 9: (Color online) Differential yields in central Pb+Pb collisions as function of transverse 
momentum pt (a,c) and transverse mass mj- (b,d) for if* (892)° (a,b) and if* (892)° (c,d) in the 
rapidity interval 0.43 < y < 1.78. The bars show statistical errors, the bands indicate systematic 
errors. Solid curves show results of exponential fits (Eq. (j3Lsee text). Dotted curves result from 
blast- wave calculations using the parameters fitted in ref. [28|. The dashed-dotted curves depict 



predictions of the UrQMD model 



29(]. 



those found for other mesons in p+p collisions 26| indicating the absence of radial flow in 
these reactions. The pr-integrated rapidity spectrum (see Fig. [12] and Table IIVI) was ex- 
tracted in the range of pt < 1-5 GeV/c, except for the last rapidity bin 1.8 < y < 2.0 which 
had a reduced range of pt < 1-2 GeV/c because of the upper momentum limit imposed 
by the dE/dx identification procedure. The range px < 1-5 GeV/c contains 99.1% of all 
K* (892) (for T = 160 MeV). The total yields (listed in Table IVTU) were extracted by fitting a 
Gaussian distribution centered at y = to the rapidity distribution. The resulting widths of 
the rapidity distributions and the mid-rapidity yields are listed in Table IIVI The extracted 
width of the rapidity distribution is consistent with the one for charged kaons 23[. The 



extracted yie 
energies (see 



ds fit well into the trend of results from p+p collisions at higher and lower 



26[ for a more detailed comparison). 
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TABLE II: Yields of K*(892)° and A"* (892)° per event in central Pb+Pb collisions as a 

function of transverse momentum pt in the rapidity interval 0.43 < y < 1.78. Both statistical 
(first) and systematic (second) errors are shown for the differential yields. The inverse slope 
parameters T of exponential fits according to Eq. (J3]) are listed with their statistical errors. 



p T [GeV/c] 


^(ir>(892) ) [(GeV/c)- 1 ] 


^(K*(892)°) [(GeV/c)" 1 ] 


0.0 - 0.5 


1.49 ± 0.12 ± 0.66 


0.64 ± 0.10 ± 0.34 


0.5 - 1.0 


3.04 ± 0.13 ± 0.38 


1.34 ± 0.09 ± 0.17 


1.0 - 1.5 


1.59 ± 0.07 ± 0.31 


0.76 ± 0.04 ± 0.17 


1.5 - 2.0 


0.55 ± 0.04 ± 0.09 


0.22 ± 0.02 ± 0.05 


T [GeV] 


0.339 ± 0.009 


0.329 ± 0.012 



TABLE III: Yields of if* (892)° and K* (892)° per event in central Pb+Pb collisions as 

a function of transverse mass tut — mo in the rapidity interval 0.43 < y < 1.78. Both statistical 
(first) and systematic (second) errors are listed. 



tut — mo [GeV/c 2 ] 


^»(^*(892)°) [(GeV/c 2 r 2 ] 


4^(^(892)°) [(GeV/c 2 )" 2 ] 


0.034 


5.95 ± 0.48 ± 2.48 


2.54 ± 0.39 ± 1.37 


0.272 


4.05 ± 0.16 ± 0.50 


1.79 ± 0.12 ± 0.22 


0.642 


1.27 ± 0.06 ± 0.24 


0.61 ± 0.04 ± 0.13 


1.070 


0.32 ± 0.02 ± 0.05 


0.13 ± 0.01 ± 0.03 



C. C+C and Si+Si collisions 

Due to the limited number of recorded events only total yields per event could be esti- 
mated (see Sect. III.C). The results with their statistical and sytematic uncertainties are 
listed in Table ELD 
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FIG. 10: Transverse momentum distribution of K*(892)° (a) and K (892)° (b) for the rapidity 
interval 0.2 < y < 0.7 in inelastic p+p collisions at 158 GeV/c incident momentum. Only statistical 
errors are shown; the overall systematic error of normalisation is 9 %. Curves show results of fits 
with the exponential function Eq. ([3]). 
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FIG. 11: Transverse mass distribution of K*(892)° (a) and K (892)° (b) for the rapidity interval 
0.2 < y < 0.7 in inelastic p+p collisions at 158 GeV/c incident momentum. Only statistical errors 
are shown; the overall systematic error of normalisation is 9 %. Curves show results of fits with 
the exponential function Eq. ([3]). 

V. DISCUSSION OF RESULTS 

Both A'*(892)° (i^*(892)°) and K + (K~) contain an anti-strange (strange) valence quark 
in addition to a light quark and should therefore be similarly sensitive to the strangeness 
content of the produced matter. The ratio of total yields (K* (892)°} /(K* (892)°} is about 2 
in C+C, Si+Si and Pb+Pb collisions (see Table IVlT]) and is similar to the ratio (K + ) / (K~) 
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FIG. 12: Rapidity distributions for K* (892)° (a) and K*(892)° (b) integrated over p T < 1.5 GeV/c 
in inelastic p+p collisions at 158 GeV/c incident momentum. Only statistical errors are shown; 
the overall systematic error of normalisation is 9 %. Curves show fits with a Gaussian function 
centered at mid-rapidity y = 0. 

~ 2.0 {2], [jij as expected. The yields per wounded nucleon are compared graphically in 
Fig. [13] (a). This quantity seems to increase from p+p to C+C and Si+Si collisions and 
then to decrease to central Pb+Pb collisions. This behaviour may result from an interplay 
between strangeness enhancement in nucleus-nucleus collisions and the interaction of the 
K* (892) and its decay products in the produced fireball. 

Because kaons and if* (892) both contain the same valence quarks, the system size de- 
pendence of the ratios (K*) / (K + ) and (K ) / (K~) is expected to be sensitive mostly to the 
interactions in the surrounding medium while the effect of strangeness enhancement should 
approximately cancel. As can be seen from Fig. [13] (b), the ratios decrease by about a factor 
2 from C+C and Si+Si reactions to central Pb+Pb collisions and about a factor of 3 when 
taking p+p reactions as the reference. Thus if* (892)° yields seem to be strongly affected by 
interactions in the produced fireball with destruction dominating regeneration. Published 
measurements from the STAR collaboration at RHIC energies 9| show a weaker suppression 
of the (K*)/ (K + ) ratio for central Cu+Cu and Au+Au collisions compared to inelastic p+p 
reactions of only about 30 %. 

Microscopic models of hadron production in nucleus-nucleus collisions have been used to 
study the modification of resonance yields during the space-time evolution of the fireball. In 
the UrQMD model ll| particle production proceeds via string excitation and decay at high 
energies, and evolves further by interactions and coalescence in the produced matter. During 
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TABLE IV: Yields of K*(892)° and K (892)° per event with statistical errors in inelastic p+p 
collisions as a function of rapidity y and integrated over px < 1.5 GeV/c. The systematic error of 
normalisation (not shown) is 9 %. Also listed are the widths a y , the central rapidity yields ^ \ v= q 
and the total yields obtained from Gaussian fits with statistical and systematic errors. 





^(ET*(892)°) 


% (K* (892)0) 


0.0 - 0.2 


n no A i 
U.0241 


± 


O.OUZO 


U.U188 


± 


U.UOzz 


0.2 - 0.4 


0.0272 


± 


0.0018 


0.0179 


± 


0.0016 


0.4 - 0.6 


0.0232 


± 


0.0013 


0.0176 


± 


0.0011 


0.6 - 0.8 


0.0192 


± 


0.0010 


0.0159 


± 


0.0008 


0.8 - 1.0 


0.0180 


± 


0.0008 


0.0146 


± 


0.0007 


1.0 - 1.2 


0.0173 


± 


0.0007 


0.0115 


± 


0.0006 


1.2 - 1.4 


0.0136 


± 


0.0007 


0.0094 


± 


0.0005 


1.4 - 1.6 


0.0113 


± 


0.0007 


0.0079 


± 


0.0005 


1.6 - 1.8 


0.0095 


± 


0.0007 


0.0048 


± 


0.0004 


1.8 - 2.0 


0.0062 


± 


0.0007 


0.0026 


± 


0.0004 


Oy 

dn 1 

dy \y=° 
total yield 


1.17 ± 0.03 ± 0.07 
0.0257 ± 0.0031 ± 0.0023 
0.0741 ± 0.0015 ± 0.0067 


1.01 ± 0.02 ± 0.06 
0.0183 ± 0.0027 ± 0.0016 
0.0523 ± 0.0010 ± 0.0047 



the model simulations track is kept of the full history for each particle thus allowing to study 
the phenomena of destruction and regeneration of resonance states 29. 30] which was first 
considered in [7J. For comparison with our measurements we extracted the predicted yields 
in the K + ir~ and K~n + decay channels, respectively, and scaled them by 3/2 to take into 
account the decay branching ratio as was done for the data. The measured rapidity spectra 
of if*(892)° and if* (892)° in central Pb+Pb collisions at 158A GeV are compared to results 
from UrQMD model calculations in Fig. [8] (dashed-dotted curves). While one observes good 
agreement for the shape of the rapidity distributions, yields are overpredicted by roughly a 
factor of two. Moreover, the predicted transverse mass distributions are steeper than those 
of the data (see dashed-dotted curves in Fig. [9] (b,d)). 

Inspection of the particle histories in the simulated UrQMD events indicates that of 
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TABLE V: Yields -J^ of K*(892)° and if* (892)° per event in inelastic p+p collisions for the 
rapidity interval 0.2 < y < 0.7 as a function of transverse momentum pr with statistical errors. 
The systematic error of normalisation (not shown) is 9 %. The inverse slope parameters T of 
exponential fits according to Eq. (|3j) are also listed with statistical and systematic errors. 



Bt [GeV/cl 


dp T dy K 1 ^ {°V*J ) 


dn (K*(892)°) 

dp T dy \^ \ ov ^) ) 


0.0 - 0.2 


0.0133 ± 0.0013 


0.0120 ± 0.0012 


0.2 - 0.4 


0.0274 ± 0.0017 


0.0227 ± 0.0014 


0.4 - 0.6 


0.0291 ± 0.0016 


0.0237 ± 0.0014 


0.6 - 0.8 


0.0171 ± 0.0012 


0.0136 ± 0.0010 


0.8 - 1.0 


0.0115 ± 0.0009 


0.00741 ± 0.00073 


1.0 - 1.2 


0.0062 ± 0.0007 


0.00466 ± 0.00050 


1.2 - 1.4 


0.0034 ± 0.0005 


0.00180 ± 0.00034 


1.4 - 1.6 


0.0016 ± 0.0008 


0.00070 ± 0.00035 


T [GeV] 


0.166 ± 0.011 ± 0.010 


0.150 ± 0.010 ± 0.010 



the originally produced AT*(892)° (A'*(892)°) about 2 % (2 %) in p+p, 12 % (13 %) in 
central C+C, 23 % (27 %) in central Si+Si and 44 % (62 %) in central Pb+Pb collisions 
are lost owing to in- medium interactions and decay (rescattering of the decay products). 
The contributions of the various mechanisms are illustrated by the curves shown in Fig. [TH 
The first reduction is due to reinteractions of the K* in the fireball medium and the second 
reduction is the effect of the decay branching ratio into K + ir~ and K~ir + , respectively. The 
last reduction is caused by the scattering of the K* decay daughters in the medium. The 
model calculations thus suggest a sizable duration of the hadronic phase of the fireball to 
allow for such reinteractions. 

Predictions of the UrQMD model for total yields are listed in Table IVIII The agreement 
with the measurements for p+p, C+C and Si+Si reactions suggests that UrQMD reproduces 
the absorption effects in these smaller systems. In contrast, predicted total A* (892)° and 
K (892)° yields for central Pb+Pb collisions exceed the experimental results by roughly a 
factor of 2. This might imply that the lifetime of the hadronic phase is larger than suggested 
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TABLE VI: Yields ^-j^^ of if* (892)° and K* (892)° per event in inelastic p+p collisions for 
the rapidity interval 0.2 < y < 0.7 as a function of my — itiq. Only statistical errors are shown. 
The systematic error of normalisation (not shown) is 9 %. 



{tut) — tuq [GeV/c 2 ] 


1 dn (K*(892)°) 
rriT dmrpdy \ \ ' ' 


1 J *\ (K*f892)°) 


0.011 


0.1331 ± 0.0132 


0.1200 ± 0.0117 


0.053 


0.0914 ± 0.0054 


0.0757 ± 0.0048 


0.132 


0.0583 ± 0.0032 


0.0477 ± 0.0027 


0.241 


0.0244 ± 0.0018 


0.0194 ± 0.0015 


0.373 


0.0129 ± 0.0010 


0.0082 ± 0.0008 


0.520 


0.0056 ± 0.0006 


0.0042 ± 0.0005 


0.680 


0.0026 ± 0.0004 


0.00138 ± 0.00026 


0.847 


0.0010 ± 0.0006 


0.00047 ± 0.00023 



by the model calculation. 

The statistical hadron gas model (HGM) was found to provide a good fit to total yields 
of stable hadrons produced in elementary e + + e~, p+p and nucleus+nucleus collisions using 
as adjustable parameters the hadronisation tem per ature T c h em , the baryo-chemical potential 
fis an d a strangeness saturation parameter 7 S 12J, |3l|. The predictions for if* (892) yields 
(which were not included in the fit of the model parameters) in p+p and nucleus+nucleus 
collisions are compared to the measurements in Fig. [15] (a). One finds that the HGM 
predictions are consistent with the measurements for p+p and light nuclei collisions, but 
exceed by more than a factor of 2 the observed yields in central Pb+Pb reactions. 

Yields of several resonance states were investigated by the NA49 collaboration, namely of 
if*(892)°, A(1520) (32! and the meson {27)]. The ratios between the measured yields and 
the predictions of the HGM model are plotted in Fig. [T3] (b) versus the respective lifetimes 
(3.91, 12.7 and 46.5 fm/c). The suppression with respect to the HGM predictions seems to 
get stronger with decreasing lifetime of the resonance. This suggests that a large part of the 
reduction of the K *(892)° yield may be caused by rescattering of its decay daughters during 
the hadronic stage of the fireball and implies that this stage lasts for a time (estimated about 
6 fm/c in Ref. 28j]) at least comparable to the lifetime of the resonance. 
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TABLE VII: Total yields of K* (892)° and K* (892)° in inelastic p+p and in central C+C, Si+Si 
and Pb+Pb collisions at 158^4 GeV beam energy. Statistical and systematic errors were added in 



29] 



)f [2]) ; 

(predictions are for the K + tt~ and K~n + decay channels, respectively, and were scaled by 3/2 to 
account for the branching ratio). 



reaction 


p+p C+C 


Si+Si 


Pb+Pb 


centrality 




min. bias 15.3% 


12.2% 


23.5% 


(N w ) 




2 14 ±2 


37 ±3 


262 ±6 


(K* (892)°^ 


this analysis 


0.0741 ± 0.0069 0.8 ±0.24 


2.2 ±0.66 10.3 ±2.5 




HGM 


0.074 0.964 


2.76 


25.1 




UrQMD 


0.076 0.74 


2.25 


22.2 


(K* (892)°^) 


this analysis 


0.0523 ± 0.0048 0.43 ±0.14 


1.3 ±0.4 


5.2 ± 1.7 




HGM 


0.041 0.455 


1.33 


12.5 




UrQMD 


0.043 0.41 


1.19 


9.5 



Alternatively, one may reconsider the assumption of simultaneous chemical freeze-out of 
the hadrons from the fireball at a unique temperature. In the statistical hadron gas model 
the yield ratio of two hadrons with identical strangeness, isospin and baryon number allows 
to estimate the freeze-out temperature Tf Q (in the Boltzmann approximation, and neglecting 
feed-down from resonances) as: 




Tfo=(m 2 - mi )/ln[ -± - 1 W 



where m 1; m 2 are the masses, Ji, J 2 the spins and A r 1 , N 2 the produced multiplicities. Insert- 
ing the numbers for the pairs A(1520), A and K*(892)°, K ± one obtains apparent freeze-out 
temperatures of 90 MeV for A(1520) and 100 MeV for K*(892)° respectively. The chemical 
freeze-out temperature fitted to the yields of stable hadrons is T c h em = 155 MeV 12J. This 
would lead to the conclusion that short-lived resonances freeze out at the end of the fireball 
evolution when the temperature has fallen below that for stable hadrons. 
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FIG. 13: (Color online) (a) Yields per wounded nucleon K* (892)° / (N w ) (dots) and 
K*(892)°/(N W ) (squares) versus size of the collision system, (b) Ratios {K*(892)°)/{K+) (dots) 
and (K (892)°)/{K-) (squares) versus size of the collision system (p+p, C+C, Si+Si and Pb+Pb 
collisions). Total kaon yields were taken from refs. 
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231 ] and appropriately scaled by {Nw)- 



For evaluating error bars the quadratic sums of statistical and systematic errors were used. 
VI. SUMMARY 

Production of the AT* (892)° and K (892)° resonances was studied via their K + tt 
and K~it + decay modes in central Pb+Pb, Si+Si, C+C and inelastic p+p collisions at 
158A GeV (t/s^ = 17.3 GeV) with the NA49 detector at the CERN SPS. Transverse 
momentum and rapidity distributions were obtained and total yields were estimated. The 
following conclusions were reached: 

• The yield of AT* (892)° exceeds that of if (892)° by about a factor of two. This ob- 
servation can be understood from the similar ratio of the K + and K~ yields and the 
valence quark composition of these mesons. 

• The yield of K*(892)° and K (892)° per wounded nucleon appears to increase from 
p+p to C+C and Si+Si collisions and then tends to decrease to central Pb+Pb reac- 
tions. This behaviour seems to reflect an interplay between strangeness enhancement 
in nucleus-nucleus collisions and attenuation of A'*(892)° and A^*(892)° in the pro- 
duced fireball. 

• The ratios (A'*(892)°)/(A" + ) and (K* (892)°) / (K~) decrease strongly with increasing 
size of the colliding nuclei. These ratios are expected to be mostly sensitive to in- 
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FIG. 14: (Color online) Predictions of the UrQMD model version 1.3pl [29| for the 23.5% most 
central Pb+Pb collisions at 158A GeV: distributions of rapidity (a,b) and transverse momentum 
p T (c,d) in the rapidity interval 0.43 < y < 1.78 of K*(892)° (a,c) and K*(892)° (b,d). The curves 
show the successive reduction of yields owing to various interaction mechanisms of the fT*(892)° 
and i^*(892)° and their decay daughters in the fireball (see text). Note that the final results show 
the yield predictions for the K + n~ and K~tt + decay channels respectively. The curves in Fig. [8] 
and Fig. [9] were obtained by scaling these by 3/2. 

teractions of the i\"*(892)° and its decay daughters with the produced dense matter. 
The decrease of the ratios suggests a substantial duration of the hadronic stage of the 
fireball. 

• The UrQMD model, although including rescattering of if* (892)° and K* (892)° and 
their decay daughters in the hadronic phase, is not able to provide a quantitative 
description of K* (892)° production in nucleus- nucleus collisions at SPS energies. 

• Yields of K* (892)° mesons in central Pb+Pb collisions are about a factor of 2.5 below 
the predictions of the statistical hadron gas model using parameters fitted to the yields 
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FIG. 15: (Color online) (a) Ratio of measured total yields of AT* (892)° (dots) and AT* (892)° 
(squares) to statistical hadron gas model (HGM) predictions [lj| versus the size of the collision 
system (p+p, C+C, Si+Si and Pb+Pb collisions), (b) Ratio of measured yield in central Pb+Pb 
collisions to the statistical hadron-gas model prediction for AT* (892)° (dot), AT* (892)° (square), cj> 
meson and preliminary measurement of A(1520) 32] versus the lifetime r of the resonance 
state. For evaluating error bars the quadratic sums of statistical and systematic errors were used. 

of stable hadrons. 

In summary, the predicted suppression of A'* (892)° yields j3] was observed in central 
Pb+Pb collisions at the SPS. It was found to be stronger at SPS than at RHIC energies. 
More comprehensive studies of the energy and system-size dependence of the suppression of 
hadron resonance production will help to better understand the hadronisation process and 
the evolution of the high-density matter droplet created in nucleus-nucleus collisions. 
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